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Abstract 

The correlated production of A and A baryons has been studied using 4.3 million multihadronic 
Z° decays recorded with the Opal detector at Lep. Di-lambda pairs were investigated in the full 
data sample and for the first time also in 2-jet and 3-jet events selected with the kx algorithm. 
The distributions of rapidity differences from correlated AA pairs exhibit short-range, local 
correlations and prove to be a sensitive tool to test models, particularly for 2-jet events. The 
Jetset model describes the data best but some extra parameter tuning is needed to improve 
agreement with the experimental results in the rates and the rapidity spectra simultaneously. 
The recently developed modification of Jetset, the Modified Popcorn Scenarium (Mops), 
and also Herwig do not give satisfactory results. This study of di-lambda production in 2- 
and 3-jet events supports the short-range compensation of quantum numbers. 
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1 Introduction 


The compensation of quantum numbers plays a key role in our understanding of the frag¬ 
mentation process whereby partons transform into observable hadrons. Consequently, baryon 
production in hadronic e + e _ annihilation final states provides data very well suited to test 
phenomenological fragmentation models. In particular, the study of di-lambda pairs allows a 
subtle testing of model predictions because of the relatively large rates and the necessity to 
compensate two quantum numbers: baryon number and strangeness. 

Fragmentation models such as Jetset jl] and Herwig [|J are based on a chainlike production 
of hadrons with local compensation of quantum numbers. In JETSET, particle production is 
implemented via string fragmentation. Baryons (B) are formed when a diquark pair is contained 
in the string (see diagram a below), thus resulting in a strong baryon-antibaryon correlation. 
This correlation can be softened by the “popcorn effect” when an additional meson (M) is 
produced between the baryon pair as shown in the diagrams b and c below. In contrast, 
Herwig describes fragmentation via the formation of clusters and their subsequent decay. 
Baryons are produced by the isotropic cluster decay into a baryon pair, which can result in 
stronger correlations than those predicted by Jetset. 
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(c) Popcorn (B(n*M)B) 
MOPS 


Di-lambda production in multihadronic Z° decays has been studied over the past years by 
experiments at Petra, Pep and Lep Q [?j|. These experiments report short-range correlations 
as observed in the distributions of the rapidities y or rapidity differences Ay of correlated 

AA pairs. The rapidity of a particle is defined as y = \ In J, where E is the energy 

of the particle and the longitudinal momentum with respect to the thrust axis. Rapidity 
differences are Lorentz-invariant under boosts along the event axis. These correlations are 
compared to predictions of Jetset and Herwig. Satisfactory agreement is found with the 
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predictions of Jetset when p, the “popcorn parameter” [], is set to the default value, p = 0.5. 
However, the tune of other parameters modeling baryon production significantly influences 
the predictions Herwig on the other hand predicts correlations much larger than those 
experimentally observed. 

The full data sample of 4.3 million hadronic Z° decays collected with the Opal detector at 
Lep in the region of the Z° peak is used in this investigation. It supplements the earlier 
OPAL work (7j by increased statistics and a more robust technique to remove the background 
contributions from the A A and AA(AA) samples in order to obtain a correlated AA sample which 
is as clean as possible. The AA correlations are investigated mainly via rapidity differences. 
They are compared to the earlier Lep results and to the predictions of Jetset and Herwig. 
The predictions of the recent Jetset modification Mops (Modified Popcorn Scenarium) |J 
are also considered. Correlated AA pairs are also studied in 2-jet events in which models can be 
tested with improved sensitivity (compared to the full data sample) when rapidity differences 
are investigated. Finally, we study correlated AA pairs within the same and within different 
jets. 

Section |2| gives a short description of the Opal detector and presents the selection of the 
A events^ in the total sample and also in 2-jet and 3-jet events, for both experimental and 
simulated data. In section |3] the separation of the AA and AA(AA) samples from the background 
and the determination of the rates of correlated AA pairs as a function of the rapidity differences 
| Ay | are discussed. Section [| contains the measured rates with their errors and a comparison to 
earlier results as well as the presentation of the differential distributions as a function of |Ay| 
and cos 9*, where 9* is the angle between the thrust axis and the A momentum calculated in the 
rest frame of the di-lambda pair. In section the models are tested using the production rates 
of A pairs as well as the cos 9* and |Ay| spectra of correlated AA pairs. The range of di-lambda 
correlations is investigated in section ^ by the assignment of the A’s to the jets. Conclusions 
are drawn in section |7|. 


2 Experimental Procedure 

2.1 The OPAL Detector 


A detailed description of the Opal detector can be found in Ref. [|T0|. Of most relevance for 


the present analysis is the tracking system and the electromagnetic calorimeter. The tracking 
system consists of a silicon microvertex detector, an inner vertex gas chamber, a large-volume 
jet chamber and specialized chambers at the outer radius of the jet chamber which improve the 
measurements in the z direction (^-chambers )[]. The tracking system covers the region | cos0| < 
0.95 and is located within a solcnoidal magnet coil with an axial field of 0.435 T. The tracking 
detectors provide momentum measurements of charged particles, and particle identification 
from measurements of the ionization energy loss, dE/dx. Electromagnetic energy is measured 
by a lead-glass calorimeter located outside the magnet coil, which covers | cos0| < 0.98. 


1 The value of p can be set in Jetset with the parameter PARJ(5): p = BB + B *l IB = 0 5+ par j( 5 ) ■ 

2 For simplicity A refers to both A and A. 

3 The coordinate system is defined so that z is the coordinate parallel to the e _ beam axis, r is the co¬ 
ordinate normal to the beam axis, <j> is the azimuthal angle around the beam axis, and 9 is the polar angle 
with respect to z. 
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2.2 Data Samples 


The analysis is based on hadronic Z° decays collected around the Z° peak from 1990 to 1995 (to¬ 
tal Lep 1 statistics). The hadronic events were selected with the standard Opal procedure [|l]J 
based on the number and quality of the measured tracks and the electromagnetic clusters and 
on the amount of visible energy in the event. In addition, events with the thrust axis close to 
the beam direction were rejected by requiring | cos thrust < 0.9, where thrust is the polar angle 
of the thrust axis. With the additional requirement that the jet chamber and the z-chambers 
were fully operational, a total of 3.895 million hadronic events remained for further analysis, 
with an efficiency of (98.4±0.4)%. The remaining background processes, such as e + e _ —> t + t~ 
and two photon events, were estimated to be at negligible level (0.1% or less). 


After the A-selection which will be described below, the selection of 2- and 3-jet events was 
performed. Charged tracks and electromagnetic clusters not associated with any track were 
grouped into jets using the kj_ recombination algorithm [I2j with a cut value y cut = 0.005. In 
addition to the standard selection criteria, the energy of the clusters and the momenta of the 
charged tracks had to be less than 60 GeV/c. To improve the quality of the jets it was finally 
required that there be at least two charged particles per jet (in addition to the possible tracks 
from A decays) and that the minimum energy per jet was 5 GeV. The cuts on the quality of 
jets were chosen to be this loose to keep the kinematic range as large as possible for comparison 
with fragmentation models. In total, samples of 1.7 million 2-jet events and 1.4 million 3-jet 
events were available for further analysis corresponding to 45% and 36%, respectively, of the 
entire data set. 


2.3 Monte Carlo Event Samples 


Monte Carlo hadronic events with a full simulation of the OPAL detector [131 and including 


initial-state photon radiation were used (a) for evaluation of detector acceptance and resolution 
and (b) for studying the efficiency of the di-lambda reconstruction as a function of the rapidity 
differences. In total, seven million simulated events were available, of which four million were 

HI 


generated by Jetset 7.4 with fragmentation parameters described in 


and three million 

tH 


The two 


were generated by Jetset 7.3 with fragmentation parameters described in 
Jetset versions differ in the particle decay tables and heavy meson resonances. There are also 
some differences in the simulation of baryon production between the two samples. Their small 
influence on the efficiency correction to the experimental data is accounted for in the systematic 
errors (see section [1.3|) . 

For comparison with the experimental results, the Monte Carlo models Jetset 7.4 and Her- 
wig 5.9 |0|] were used. Both models give a good description of global event shapes and many 
inclusive particle production rates, but differ in their description of the perturbative phase and 
their implementation of the hadronization mechanism. 


Tracks and clusters are selected in the Monte Carlo events, which include detector simulation, in 
the same way as for the data, and the resulting four-vectors of particles are referred to as being 
at the ‘detector level’. Alternatively, for testing the model predictions, Monte Carlo samples 
without initial-state photon radiation nor detector simulation are used, with all charged and 

4 The fragmentation parameters of Herwig 5.9 were identical to those used in our tuned version of Her- 
WIG 5.8 @ with the exception of the maximum cluster mass (CLMAX) which was set to 3.75 GeV in order to 
improve the description of the mean charged particle multiplicity in inclusive hadronic Z° decays. 
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neutral particles with mean lifetimes greater than 3x10 10 s treated as stable. The four-vectors 
of the resulting particles are referred to as being at ‘generator level’. 


2.4 A Reconstruction 


m 


7T p as described 

Briefly, tracks of opposite charge were paired and regarded as a secondary vertex 


Neutral strange A baryons were reconstructed in their decay channel A 

m 


candidate if the track pair intersection in the plane perpendicular to the beam axis satisfied 
the criteria of a neutral two-body decay with a decay length of at least 1 cm. 


Each candidate track pair was refitted with the constraint that the tracks originated from a 
common vertex, and background from photon conversions was suppressed. Information from 
dE/dx measurements was used as in [[L7[] to help identify the 7r and p for further background 
suppression, primarily due to Kg^ 7r + 7r - . Two sets of cuts, called ‘method 1’ and ‘method 2’ are 
described in JT7J for A identification. For the present analysis, A candidates were reconstructed 
using method 1, which is optimized to have good mass and momentum resolution. 


By these means a narrow A mass peak above a small background has been ob¬ 
tained. The selection of di-lambda candidates with both invariant masses in the range 
1.1057 GeV/c 2 < m^p < 1.1257 GeV/c 2 (region A in figure jlj) retains most of the A signal 
for further analysis. 


3 Selection of Correlated A-pairs 

3.1 Method 

Events with more than one A candidate that had passed the above selection criteria were 
considered and all possible pair combinations of the A and A baryons within an event were 
formed. This resulted in pairs of AA, AA and AA. Combinations were rejected if the pair had 
a track in common. The remaining pairs are henceforth referred to as A-pair candidates. 

The three types of baryon pairs can be grouped into two classes: pairs with different baryon 
numbers A A and pairs with equal baryon numbers AA(AA). Only in A A pairs can the baryon 
and flavor quantum numbers be compensated by correlated production. AA(AA) pairs can 
never be produced in correlation and hence they will occur only in events with more than one 
baryon-antibaryon pair (BB). In such events uncorrelated AA pairs from different (BB) pairs 
are also possible. The number of uncorrelated AA pairs corresponds to the number of pairs 
with same baryon number. Hence, the number of correlated AA pairs can be derived via 

KT=n aa -(n aa + n aa ). (i) 

At this stage 9479 AA and 4217 (AA+AA) pair candidates are selected. 

3.2 Background Subtraction and Efficiency Correction 

Due to the small statistical errors it is necessary to keep systematic uncertainties as low as pos¬ 
sible in this analysis. The correct subtraction of non-A background from the pairs is therefore 
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of particular importance. This background consists mainly of other long-lived particles with 
similar decay topologies (namely Kg—> 7r + 7r _ ) and random track combinations. An important 
contribution to the contamination is the so-called correlated background from Acandidates that 
have been reconstructed with one false decay track. They are more numerous in pairs with op¬ 
posite baryon number because the number of AA pairs is far higher than the number of AA(AA) 
pairs. For this reason the background has to be estimated in the two samples separately. Back¬ 
ground pairs occur when either one or both A-canclidates are fake. In the two-dimensional mass 
plane in figure [I], pairs with one fake A form horizontal and vertical bands of background, while 
pairs with two fake candidates are uniformly distributed in the region above the lower mass 
bounds. 

The background was subtracted using a two-dimensional sideband method. The background in 
the signal region A was measured from two mass windows (sidebands) of the same size (regions 
Bj and B 2 ) placed in the two bands of background. In this way the background with two 
fake candidates is counted twice. The latter was determined from region C outside the bands. 
Hence, the signal is obtained with the subtraction: 

Signal = N a - (N Bi + N B2 - N c ) . 

We optimized the position of the sidebands with a MC test investigating the deviations between 
the background-corrected sample and the true-A sample. The stability of this method was tested 
in the experimental data by shifting the position of the sidebands by one half of the band size 
from the optimized position. The fluctuations were of the same size as the deviations found in 
the MC. 

Finally the background-corrected A A and AA(AA) signal distributions were corrected for de¬ 
tector acceptance and reconstruction efficiency as functions of Ay and cos#*. The average 
efficiency in the total hadronic sample was found to be ~ 2%, varying between 1.3% and 2.5% 
over the |Ay|/cos#* range. 


4 Experimental Results 

In an earlier OPAL paper [[7] based on the 1990 and 1991 data samples we already investigated 
the production dynamics of baryon-antibaryon pairs. In this section, we present the rates and 
differential distributions of A pairs using the full 1990 to 1995 LEP 1 data in three samples: 
the entire set of multihadronic events, the 2-jet and the 3-jet events. 

4.1 Pair Production Rates 

The resulting rates for A A and AA(AA) pairs in all hadronic events, determined as sum over 
all corrected |Ay| bins, are given in table [I]. The rates for the correlated A A pairs are derived 
according to equation (1) from the difference of the opposite and same baryon number pairs. 
Compared to the results from other Lep experiments and to the previous OPAL publication, 
good agreement is found. 

The di-lambda rates in 2- and 3-jet events are listed in table ||. In 3-jet events, due to the 
higher color charge of the gluons, the average pair multiplicity is higher. 


4.2 Differential Distributions 


We studied the correlations in the differential AA spectra using the observables |Ay| and cos 9* 
as they are particularly sensitive for comparison with Monte Carlo models. The differential 
distributions are shown in figure |2j. The short-range correlations show up as a peak in the 
region \Ay\ < 2.0. 

When investigating |Ay| distributions, we will restrict ourselves to 2-jet events. This is due to 
the fact that in 3-jet events many particle momenta have large angles to the thrust axis, resulting 
in smaller longitudinal momenta and smaller rapidity differences, independent of correlations. 
As a result the Ay distribution is broader and less steep in 2-jet events than in the 3-jet or 
the total sample (see figure ||a). Consequently, also the range of variations is larger in 2-jet 
events and yields a higher sensitivity in the comparison with model predictions. 


4.3 Systematic Errors 

The systematic error is found to be largely independent of |Ay| and cos#*, and in the subse¬ 
quent discussion of the differential distributions of the correlated pairs, only normalized dis¬ 
tributions are considered. These are largely insensitive to effects of systematic uncertainties. 
Consequently, the systmatic errors discussed below are mainly relevant for the total rates. 

For the determination of the experimental uncertainties we considered the following sources of 
systematic effects: 


• Uncertainties due to the subtraction of background via the sidebands. These were es¬ 
timated using simulated events by applying the analysis to the fully detector simulated 
MC and comparing the rate from the background corrected sample to the true number. 

• Efficiency uncertainties. These were estimated from the difference of the results when 
the efficiency correction was done using both Jetset versions 7.3 and 7.4 samples in 
combination and using them separately. 

• The statistical error of the efficiency due to the limited sample size of the simulated events 
at detector level. 


Uncertainties in the modelling of the cut variables used for the A selection. This error is 
taken from a former analysis Q where it was determined very precisely for single A’s. 
The error given there is doubled for the A pairs in the present analysis. 


These effects contribute to the total systematic error as shown in table ^], where the relative 
systematic errors from the different sources are compared to the total systematic as well as to 
the statistical error. Statistical and total systematic errors contribute about equally. 


5 Comparison with Fragmentation Models 

We start the discussion with the numbers and distributions of the models with OPAL default 
tunes that optimize the general performance of the models and the agreement with the measured 
single particle rates. 
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5.1 Pair Production Rates 


We investigate the di-lambda rates first in the total hadronic data sample comparing the mea¬ 
sured rates to the predictions of the models Jetset 7.4, Mops and Herwig 5.9 (see table |l|). 
None of the models gives a perfect description of the data but HERWIG clearly exhibits the 
largest disagreement. 


The comparison of the di-lambda rates in 2- and 3-jet events is given in table |2j. The higher 
multiplicity in 3-jet events compared to 2-jet events is qualitatively well described by all three 
models. However, only JETSET yields a prediction compatible with the measured numbers. In 
the 2-jet event sample the agreement is excellent. In the 3-jet sample all the measured rates 
exceed the Jetset predictions. This can be compared to the observation that A rates in gluon 
jets are too low in Jetset JT8[ . 


5.2 Differential Distributions 

To further investigate the nature of the AA correlations we compare the differential distributions 
of correlated AA pairs with the predictions of the various models. We use the variables cos 9* 
and | Ay | and test their sensitivity to distinguish between the different fragmentation models 
and baryon production mechanisms. All distributions are of the type d ^ y \y N being the 
total number of entries. This has the advantage that they are independent of the total rates 
and that the systematic errors mostly cancel out, since they are nearly independent of both 
| Ay | and cos 9*. 

The angle 9* is particularly suited to distinguish between string and cluster fragmentation. The 
mostly isotropic cluster decay (Herwig) results in a relatively flat cos 9* distribution whereas 
string fragmentation produces the correlated AA system predominantly close to the thrust axis, 
i.e., with cos 9* rs 1. These predictions are compared to the measurement in figure The data 
show a distribution that is strongly peaked towards 0056 **= 1 and therefore clearly rule out the 
Herwig cluster model. The predictions of MOPS agree somewhat better with the experimental 
distribution but they also fail to model the forward peak correctly. Only Jetset yields a good 
description of the data. 

On the other hand, especially in 2-jet events, the rapidity difference |Ay| is more sensitive 
to show differences in the strength of the correlations. The experimental data and model 
predictions are compared in figure (|. Again Jetset gives the best, albeit not completely 
satisfactory, description of the measured distribution. Herwig generates correlations which 
are far too strong . The MOPS model with its built-in facility to allow for several “popcorn 
mesons” should yield weaker correlations than Jetset; however, in contrast to this naive 
expectation it produces a narrower |Ay| distribution, i.e. stronger correlations. We see the 
following possible reasons for this: first of all, and different from JETSET, a new kinematic 
property is built into MOPS: the low-T-suppression 0. This suppresses popcorn fluctuations 
at early times in the color held, resulting in very strong correlations. Secondly, it appears that 
the strength of the correlations is influenced more by the rate of baryon production via the 
popcorn mechanism than by the actual number of intermediate mesons produced. As JETSET 
and MOPS are tuned to show the same mean number of popcorn mesons instead of popcorn 
systems, MOPS has fewer popcorn systems and therefore stronger correlations. 
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5.3 Tuning of Models 


In an earlier Opal analysis of strange baryons [|] it was observed that the agreement be¬ 
tween experimental data and Jetset predictions can be improved by adjusting some of the 
diquark parameters: improving the predicted shape of the |Ay| distribution was possible by 
varying the popcorn parameter, p=PARJ(5), that influences the frequency of popcorn produc¬ 
tion and hence the correlation strength. It was found that p acts on both the shape of the 
rapidity spectrum and the production rates. Two other parameters were used to correct for 
this change of predicted multiplicities: the ratio of the strange to non-strange diquarks over 
strange to non-strange quarks, (us:ud/s:d) = PARJ(3), and the ratio of spin-1 to spin-0 di¬ 
quarks, (1/3 ■ [qq]i/[qq] 0 )=PARJ(4). These last parameters affect mainly the rates and leave 
the spectra nearly unmodified. When attempting to improve the predictions of the MOPS model 
in the same manner, the most direct correspondence to p in JETSET is the MOPS parameter 
PARJ(8)=/3(u), the transverse mass of an intermediate u-quark. The higher the transverse 
mass of the intermediate system (with several quark pairs possible), the lower the probability 
to produce this popcorn system and the stronger the correlations. Therefore, in both JETSET 
and MOPS we tried first to improve the agreement with the data distributions by tuning the 
parameters that influence the correlation strength (figure |5| for JETSET.) In Jetset, the pop¬ 
corn probability p was varied from 0%-90%, while in MOPS the transverse mass of a u-quark, 
/3(u)=PARJ(8), was altered between 0.2 and 1.0 GeVW All other parameters remained at the 
Opal default values. As can be seen from table ||, for the p parameter, these variations affect 
not only the shape of the |Ay| distribution but also the di-lambda rates, as expected. 

The predictions with the different popcorn parameter values in Jetset are compared to the 
data in figure ^a. Only the results from parameter settings above the default value of p = 0.5 
are shown, since lower values give a poorer agreement with the data. The best agreement is 
found in the range 0.6 < p < 0.8. Popcorn values within this range also yield good agreement 
between data and predictions for the cos 6* distribution. However, when the influence of the 
popcorn parameter on the predicted di-lambda rates is also considered (table pj) use has to be 
made of the other two JETSET parameters that affect the strange baryon production in order 
to tune the rates back to values corresponding to the measurement. It can be seen from table [5] 
and figure [5|b that such a tune clearly produces a better agreement with the rates (also for single 
particle production) while it does not change the spectra of rapidity differences significantly. 
Using the results of other Opal analyses, it can also be seen that the tune does not change 
the strange meson (Kg) rate, nor does it affect the non-strange baryon (p) rate significantly. 
The known problems [|3] in modeling the decuplet baryon rates are also seen here. No further 
attempt has been made to globally optimize the parameter set, however. 

The tune of parameter PARJ(8) in MOPS did not result in an improvement. Although the value 
of the parameter was varied in a comparatively wide range, the effect on the |Ay| distribution 
was almost imperceptible. PARJ(8) clearly is not suited to adjust the Mops model to the data. 
Therefore, we tested another parameter of the model using the relative difference between the 
fragmentation function f(z ) for baryons and mesons, the parameter PARJ(45). Again, the 
variation did not notably change the shape of the distribution. This relatively poor perfor¬ 
mance of the Mops Monte Carlo in describing the Ay dependent AA correlations seems to be 
connected to the known shortcomings of the model in describing pj_-related distributions [pj|. 
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6 Di-lambdas in Jets 


After studying the strength of AA correlations in the Ay spectra, we will now present results 
on the range of the correlations by assigning both partners from a correlated pair to the recon¬ 
structed jets in an event. For short-range correlations both partners are expected within the 
same jet whereas long-range correlations (which can be obtained by the production of baryons 
from the primary quarks) should result in an assignment to different jets. We use the follow¬ 
ing two classifications for the assignment study: both partners within the same jet, and each 
partner in a different jet. 

Due to the fact that it is impossible to map 2-jet events at detector level to 2-jet events at 
generator level, we do not attempt to apply efficiency corrections but compare our uncor¬ 
rected results with the Jetset predictions at detector level. We count the number of AA and 
AA(AA) pairs in each sample and obtain the number of correlated pairs again from the rela¬ 
tion _/\T^ relatcd = A aa — (Aaa + Ayy). The amount of background in like- and unlike-sign pairs 
approximately cancels out in this subtraction as long as the contribution from the correlated 
background (see section 3.2) can be neglected. The numbers of pairs obtained from the same 
jet and from different jets are listed in table [4| for both 2- and 3-jet events. The major part 
of the correlated pairs is reconstructed within the same jet (about 96% in 2-jet events, 81% 
in 3-jet events) whereas only a very small fraction is found in different jets. These experi¬ 
mental numbers are in excellent agreement with the Jetset predictions at detector level and 
support the assumption of short-range compensation of baryon number and strangeness in the 
fragmentation process. 


7 Summary 

AA correlations have been studied in 4.3 million multihadronic Z° decays, with the correlated 
sample obtained from the difference: AA corr =AA-(AA+AA). The analysis has been performed 
in terms of cos 6* and rapidity differences |Ay|. As the rapidity is defined with respect to 
the event (thrust) axis, the sensitivity of the analysis is seen to be higher in 2-jet events. 
Therefore three data samples have been analyzed: the entire hadronic event sample, 2-jet 
events (45%), and 3-jet events (36%). The experimental findings have been used to study the 
baryon production mechanism implemented in various phenomenological fragmentation models. 

The following results have been obtained: 

• In the full data set, the measured production rates of AA, AA(AA) and, consequently, 
AA corr are in good agreement with a previous OPAL measurement and results from Aleph 
and DELPHI, but show significantly smaller errors. 

• The cos 6* distribution of correlated A-pairs is well suited to distinguish between isotropic 
cluster and non-isotropic string decay, and clearly favours the latter, implemented in 
JETSET. The predictions of the isotropic cluster model HERWIG are ruled out by the 
data: they do not describe the features of correlated AA production. 

• The rapidity difference |Ay| is used to study the strength of correlated di-lambda pro¬ 
duction. The measured distribution exhibits strong local correlations. 
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• Satisfactory reproduction of the experimental results is obtained with the predictions of 
the string fragmentation model JETSET. Improved agreement can be found by tuning 
some of the default parameters used by Opal. After adjusting the popcorn parameter, 
to improve the description of the |Ay| spectrum, other parameters, fixing the fraction of 
diquarks with strangeness and spinl, have to be modified to readjust the predicted rates 
to the experimental ones. This procedure does not affect the previously optimized Ay 
distribution. The Herwig model cannot describe the measured |Ay| spectra, and the 
predictions of MOPS, a recently published modification of the JETSET model, also fail to 
reproduce the experimental data, even after some parameter tuning. 

• In the 2-jet and 3-jet event samples it is found that correlated AA pairs are produced 
predominantly within the same jet, supporting the assumption of a short-range compen¬ 
sation of quantum numbers. Again, the JETSET predictions are in good agreement with 
the experimental results. 

In conclusion, the analysis of correlated di-lambda pairs proves to be a very effective tool to 
test fragmentation models. Jetset is the only candidate model studied, describing the data 
successfully. 
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Tables 



Apairs/Hadronic Event [ xlO 2 ] 

AA 

AA(AA) 

-^-A corr 

This Analysis 

8.95 ± 0.15 ± 0.31 

2.83 ± 0.11 ± 0.17 

6.12 ±0.19 ± 0.28 

Opal J7[ 


8.26 ± 0.42 ± 0.79 

2.05 ± 0.39 ± 0.28 

6.21 ±0.54 ± 0.84 

Aleph | 


9.3 ± 0.9 

2.8 ± 0.3 

6.5 ± 1.0 

Delphi 

| 

9.0 ± 0.9 

1.8 ± 0.6 

7.2 ± 1.1 

Jetset 7.4 

7.75 

2.24 

5.51 

Mops 

10.57 

2.63 

7.94 

Herwig 5.9 

15.09 

3.06 

12.03 


Tabic 1: Comparison of average A pair multiplicities from this analysis with those from a pre¬ 
vious OPAL analysis, with the results from other Lep experiments and with model predictions. 
The statistical error is given first, the systematic error second. For Aleph and DELPHI only 
the total error is available. 



AA 

AA(AA) 

AA CO rr 


Ap a i rs /2-Jet Event [ xlO 2 


Opal data 

5.99 ± 0.21 ± 0.30 

1.44 ± 0.14 ± 0.15 

4.55 ± 0.25 ± 0.31 

Jetset 7.4 

6.14 

1.45 

4.69 

Mops 

8.34 

1.68 

6.66 

Herwig 5.9 

13.16 

2.45 

10.71 


Apairs/3- Jet Event xlO 2 


Opal data 

9.55 ± 0.24 ± 0.41 

2.98 ± 0.18 ± 0.22 

6.67 ± 0.30 ± 0.32 

Jetset 7.4 

8.70 

2.66 

6.04 

Mops 

11.92 

3.17 

8.75 

Herwig 5.9 

16.11 

3.31 

12.80 


Table 2: Average multiplicity of A pairs in 2- and 3-jet events compared to model predictions. 
The statistical error is given first, the systematic error second. 
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Source of Error 

Effect on the AA corr Rate 

All Hadr. 

2-Jets 

3-Jets 

Background Systematics 

0.9% 

2 .8% 

0.9% 

Jetset 7.3/7.4 Mixing 

2.3% 

3.5% 

1.1% 

Monte Carlo Statistics 

2.3% 

4.3% 

3.5% 

Cut Simulation 

3.0% 

3.0% 

3.0% 

Total Syst. Error 

4.5% 

6 .8% 

4.8% 

Stat. Error 

3.1% 

5.5% 

4.5% 


Table 3: Relative errors in measuring the multiplicity of correlated A pairs in the three event 
samples. 


Assignment 

AA 

AA(AA) 

-A-A-corr 

Fractio 

Opal data 

n of AA corr 

Jetset det. level 

2-Jet Events 

Same 

Jet 

1994 

469 

1525 

(95.6 ± 2.3)% 

(95.8 ± 1.8)% 

Different 

Jets 

719 

649 

70 

(4.4 ± 2.3)% 

(4.2 ± 1.8)% 

3-Jet Events 

Same 

Jet 

2088 

409 

1679 

(80.6 ± 1.8)% 

(80.6 ± 1.4)% 

Different 

Jets 

1174 

769 

405 

(19.4 ± 1.8)% 

(19.4 ± 1.4)% 


Table 4: Assignment of A pairs to the reconstructed jets in 2-jet and 3-jet events, compared to 
the predictions of Jetset 7.4. The errors are statistical only. 
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Jetset 7.4 

Data 

p 

0.5* 

0.7 

0.7 

PARJ(3) 

0.45* 

0.45* 

0.60 

PAR J (4) 

0.025* 

0.025* 

0.010 

Di-lambda Pairs in the Total Hadronic Sample 

AA 

0.0775 

0.0668 

0.0859 

0.0895 ± 0.0034 

AA(AA) 

0.0224 

0.0187 

0.0256 

0.0283 ± 0.0020 

-A-A-corr 

0.0551 

0.0481 

0.0603 

0.0612 ± 0.0034 

Di-lambda Pairs in 2-Jet Events 

AA 

0.0614 

0.0518 

0.0677 

0.0599 ± 0.0037 

AA(AA) 

0.0145 

0.0119 

0.0164 

0.0144 ± 0.0021 

-A--A-cor r 

0.0469 

0.0399 

0.0513 

0.0455 ± 0.0040 

Particle 

Multiplicities 

K° 

2.02 

2.03 

2.02 

1.99 ± 0.04 

proton 

0.93 

0.89 

0.92 

0.92 ± 0.11 

A 

0.338 

0.316 

0.361 

0.374 ± 0.010 

£+ 

0.075 

0.067 

0.087 

0.099 ± 0.015 

£° 

0.073 

0.065 

0.086 

0.071 ± 0.018 

£- 

0.068 

0.059 

0.080 

0.083 ± 0.011 

— 

0.0278 

0.0241 

0.0341 

0.0259 ± 0.0011 

A++ 

0.10 

0.12 

0.08 

0.22 ± 0.06 

£(1385)± 

0.0457 

0.0546 

0.0393 

0.0479 ± 0.0044 

“(1530)° 

0.0036 

0.0040 

0.0035 

0.0068 ± 0.0007 

o- 

0.0006 

0.0006 

0.0006 

0.0018 ± 0.0004 


Tabic 5: Comparison of inclusive di-lambda yields with Jetset Monte Carlo predictions using 
the Opal default tune (second column), a tune to obtain agreement with the measured |Ay| 
spectra (third column) and a tune to obtain simultaneous agreement in distributions and rates 
(fourth column). The measured values with total errors are given in the fifth column. The single 
particle inclusive rates are given for comparison: the experimental numbers for Kg, protons, 

0. respectively; the remaining numbers 


19 , 20 , 21 


£ and A ++ baryons are taken from 

are from [|H]]. The parameters used for the tune are described in the text. Parameter default 
values are marked with a star. 
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Figures 



Figure 1: Two-dimensional mass distribution of AA candidates and projections onto the mass 
axes. The background forms a horizontal and a vertical band from pairs with one fake A 
above a uniform background from two non-A candidates. The signal peak at the A mass of 
1.116 GeV/c 2 is clearly visible. 
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1/N dN/d(cos 0*) 1/N dN/d(IAyl) 




cos 0* 


Figure 2: Comparison of the shape of differential distributions in all hadronic events and in 2-/ 
3-jet events for the Ay distribution in (a) and the cos 9* distribution in (b). The errors shown 
are purely statistical, the influence of the systematic errors is negligible. 
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Figure 3: Comparison of the measured distribution of the angle 6* for correlated A A pairs in 
all hadronic events with the predictions from the various models. The errors shown are purely 
statistical, the influence of the systematic errors is negligible. 
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lAyl 


Figure 4: Comparison of the rapidity difference distribution of correlated AA pairs from the 
2-jet events with the model predictions.The errors shown are purely statistical, the influence of 
the systematic errors is negligible. 
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lAyl 


Figure 5: Comparison of the measured rapidity difference distribution of correlated A A pairs in 
2-jet events with the JETSET predictions, (a) Different values of the popcorn parameter (while 
the other parameters remain at their default values) (b) The best popcorn value (0.7) with and 
without a retune of PARJ(3,4). The errors shown are purely statistical, the influence of the 
systematic errors is negligible. 
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